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The mean free paths and in-medium scattering cross sections of energetic nucleons in neutron-rich 
nucleonic matter are investigated using the nucleon optical potential obtained within the relativistic 
impulse approximation with the empirical nucleon-nucleon scattering amplitudes and the nuclear 
densities obtained in the relativistic mean field model, ft is found that the isospin-splitting of 
nucleon mean free paths, sensitive to the imaginary part of the symmetry potential, changes its sign 
at certain high kinetic energy. The in-medium nucleon-nucleon cross sections are analytically and 
numerically demonstrated to be essentially independent of the isospin asymmetry of the medium and 
increase linearly with density in the high energy region where the relativistic impulse approximation 
is applicable. 
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I. INTRODUCTION 

The in-medium nucleon-nucleon (NN) cross sections and the associated nucleon mean free paths 
(MFP) in neutron-rich nucleonic matter are important for understanding both the structure of rare 
isotopes and the reaction dynamics involving these nuclei. For instance, much information about 
the size of neutron skins and nucleon density profiles of neutron-rich nuclei can be obtained from 
Glauber model analyses of the total reaction cross sections[l]. A key input in these analyses is the 
isospin dependent in-medium NN cross sections. Usually, however, in most Glauber calculations 
only the isospin averaged free-space nucleon-nucleon cross section is used. The knowledge about 
the in-medium NN cross sections in neutron-rich matter will thus affect significantly the extracted 
structure information about rare isotopes. It is also critical for understanding nuclear reactions as 
the in-medium NN cross sections determine the stopping power and generally transport properties of 
nuclear matter. The in-medium NN cross sections, as basic inputs of transport models, are important 
in determining essentially all dynamical observables, such as, the collective flows, isospin diffusion, 
energy and momentum transfers, and particle production [2-7]. While the in-medium NN cross 
sections and the nuclear equation of state (EOS) are in principle determined by the same nuclear 
effective interaction, although they are sensitive to different parts and ranges of the interaction in the 
medium at different densities, much more attention has been putting on investigating the EOS than 
on the in-medium NN cross sections in neutron-rich matter. However, it was demonstrated recently 
that the major theoretical uncertainty in extracting the density dependence of the symmetry energy 
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of neutron-rich matter from heavy-ion reactions is our poor knowledge about the isospin-dependent 
in-medium NN cross sections [3]. It is encouraging to notice that more work is being done recently in 
calculating the in-medium NN cross sections using various approaches, such as, the relativistic and 
non-relativistic Brueckner theory [5, 8-16] and the Green function approach [17-19]. Nevertheless, 
most of these calculations are for in- medium NN cross sections in symmetric nuclear matter. Very 
little information is available for the isospin dependence of the NN cross section in neutron-rich 
matter. While the available calculations are indeed very useful, necessarily are the NN cross sections 
and the EOS consistently calculated using the same nuclear effective interactions. 

The nucleon-nucleus optical potential plays a central role in determining the in-medium NN cross 
sections and the nucleon MFP. In particular, the imaginary part of the optical potential allows 
one to extract easily the MFP and thus also the in-medium NN cross sections. There are sev- 
eral possible ways to derive the optical potential [20]. It can be calculated from various mod- 
els either non-relativistically or relativistically, with the microscopic Brueckner approach, see, e.g., 
Ref.[10, 12, 13, 21, 22]. Theoretically, approximations have to be made to render the calculations 
feasible. For instance, the evaluation of the imaginary part of the optical potential depends on the 
treatment of the widths of intermediate states that are rather unclear or the nucleon polarization in 
the medium that is in principle coupled to the vacuum. On the other hand, one can start from a 
physically reasonable approximation for the optical potential and determine its parameters using the 
experimental data [23]. The relativistic impulse approximation (RIA) actually follows this way. In 
the RIA, the optical potential is obtained in a form similar to the non-relativistic tp approximation. 
Within the RIA, the basic ingredients of the optical potential are the free Lorentz invariant NN scat- 
tering amplitudes and the nuclear scalar and vector densities in nuclear matter [24-28] . An attractive 
feature of the RIA is that the relativistic optical potentials are experimentally constrained by the 
free-space NN scattering data. The nuclear densities are calculated from the relativistic mean-field 
(RMF) models that provide also a dynamical description for the spin-orbit coupling [29]. Along with 
the success of the RMF models in describing the nuclear structure, the RIA was justified by nicely 
reproducing the proton-nucleus elastic scattering data at high energies. Of course, the limitation 
of the RIA is that it is valid only at reasonably high energies in not so dense matter. As most of 
the existing microscopic calculations were devoted to the low and intermediate energy regions, the 
in-medium NN cross sections at high energies remain largely unknown. Our study here of the nucleon 
MFP and in-medium NN cross sections at high energies within the RIA is thus complementary to 
the existing work. Moreover, given the exciting scientific opportunities provided by the constructions 
of several high energy radioactive beam facilities around the world, the knowledge on the mean free 
paths and in-medium scattering cross sections of energetic nucleons in neutron-rich nucleonic matter 



is certainly necessary and useful. 

The paper is organized as follows. In section II, the main formulas of the RMF, MA, and the 
calculations of the nucleon MFP's and in-medium NN cross sections are briefly reviewed. The results 
and discussions are presented in section III. Finally, a summary is given in Section IV. 

II. A BRIEF SUMMARY OF THE FORMALISM 

In the present work, we use the relativistic model Lagrangian [30] 

C = #7^ — M + g*a - g^uf - g* pl ^ - e±(l + r 3 )^A^ + ±(8^8"* - mfa 2 ) 

—V" + \B, V B^ + \mfb X \A^\ (1) 

where ip,a 7 uj, and b are the fields of the nucleon, scalar, vector, and isovector-vector mesons, with 
their masses M, m*,m*, and m*, respectively. The meson coupling constants and masses with 
asterisks denote the density dependence, and they are given by the BR scaling [30-32]. The Dirac 
spinor in momentum space is written as 



, E* + M* ( 1 \ 

with M* = M - g*a, E* = y/M* +k 2 , and X s the Pauli spinor. The (Oth-component) vector and 
scalar densities are respectively give by 



with i = p,n. With Dirac spinors, the NN scattering amplitudes can be defined [28]. The Lorentz 
invariant NN scattering amplitude can be decomposed into five components 

T = F s + F V 7f 72„ + F T a? v a 2flv + Fp^l + i^TiT^I^. (5) 

where S,V,T,P and A denote the scalar, vector, tensor, pseudoscalar and axial vector, respectively. 
The five components are dependent on the squared momentum transfer q 2 . The scattering amplitude 
in the RIA is obtained from the expectation value of a sum of single two-body NN scatterings. For the 
spin-saturated nucleus, only survived are the scalar (Fs), vector (Fy) and tensor (Ft) terms, while 
the scalar and vector parts dominate the scattering amplitude. In the RIA, the optical potential in 
momentum space is thus the sum of the products of the density and amplitude of respective scalar 
and vector components [24, 25] 

UoM = -^^[F s (q)/5s(q) + 70*V(q)p B (q)], (6) 



where pi a b and M are the laboratory momentum and mass of the incident nucleon, respectively. This 
expression is similar to the tp approximation in the non-relativistic approximation. The momentum 
dependence of the optical potential can give details of the angular distribution of nucleon scatterings, 
while the total cross section depends only on the forward scattering where q = 0. Following Rcf. [33], 
here we are interested in the optical potential in nuclear matter where the densities in coordinate 
space are constant. This means that a delta function <5(q) is induced so that q = in the scattering 
amplitudes in Eq.(5). Therefore the optical potential in nuclear matter in coordinate space is given 
as: 

U opt = -^±[F s p s + l0 F v p B }, (7) 

where ps and ps are the spatial scalar and vector densities of nucleons, respectively. For the 
application to the scattering on finite nuclei, one can resort to the local density approximation. 
Similar to Ref. [33], here we adopt the F$ and Fy determined directly from the experimental NN 
phase shifts [28]. The optical potential is decomposed into the following form 

Uopt^U^+joU^. (8) 

with Ug ot and J7g°* denoting, respectively, the scalar and (Oth-component) vector parts of the optical 
potential. Since Ug ot and J7q°* are generally complex, they can be further expressed as 

V% ot = U s + iW s , Ul ot = U + iW . (9) 

An application of rclativistic optical potentials is to evaluate the nucleon MFP based on the 
dispersion relation which in the relativistic frame is written as 

{E k -Ul ot f = k 2 + (M + U t s ot ) 2 . (10) 

Equivalcntly, it can be written in a form similar to the non-relativistic approximation [34] 

k 2 k 2 

2M = 2M +U ^ {Ekm) > (U) 
where the Schrodinger equivalent potential (SEP) is given by 

1 Tjtot 



U\% = Ur + Ul ot + ^(Ur 2 - Ul ot2 ) + -jfE kin , ( 1 2 ; 



with fc^ = El m + 2ME km , E k = E kin + M, and = U sep +iW sep . The real part of the symmetry 
potential is 



with the isospin asymmetry 5 = (p n — p p )/pB- The above U sym taken at the normal density is the 
well known Lane potential extracted from nucleon- nucleus scatterings [35]. Similarly, we can define 
an imaginary symmetry potential as 
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W aym = sep ns sep . (14) 



Using the dispersion relation (11) and assuming a complex momentum k = k^ + iki, the nucleon 
MFP A is derived exactly as [36] 

A = = \[-M{E Un + Hp - U sep ) + M({E km + fff - U sep ) 2 + Wi p y/T 1/2 - (15) 

A much simpler expression for the nucleon MFP can be derived by expanding the momentum in the 
vicinity of kn [37]. The real and imaginary parts of the momentum is thus obtained as 

k R = (E 2 kin + 2ME kin - 2MU sep ) 1 '\ k T = -W sep (^ + (16) 

Since there is no explicit momentum dependence in the optical potentials with the RIA, the nucleon 
MFP is written as 

A ^ip~^fe;' z=p ' n - (17) 

Numerically, we noted that almost the same nucleon MFP can be obtained from either Eq.(15) or 
(17), while the latter is more convenient for making analytical analyses. The nucleon MFP can also 
be measured as the length of the unit volume defined by the matter density and the NN cross section, 
and it is thus expressed by [38] 

Ai = {p p (T* p + p„CT*„) _1 , (18) 

where p p and p n are respectively the proton and neutron densities. The in-medium NN cross sections 
can be obtained by inverting the above equation. To write the results compactly, we may define the 
following two quantities, 

2 \ n Ap 25 \ n Ap 

They can be further written in terms of the imaginary parts of the symmetry potential and the 
isoscalar SEP using Eqs.(14) and (17) 

A_1 = l^ Wsep ' ~ X l = T^ Wsymi (20) 

where the k r ^ p are approximated by the fc#. This is a very good approximation because at high 
energies the U sym is negligible compared to the kinetic energy and the isoscalar SEP's given by 

Use P = (U? ep + U! ep )/2, W sep = {Wf ep + Wf ep )/2. (21) 
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This point will be shown numerically in the next section. As it was pointed out in Ref. [33], the 
U sym itself is isospin independent because the difference between the neutron and proton potentials 
is largely linear in isospin asymmetry. We find that such an isospin independence retains also 
for the Wsym- It is therefore understandable that the A and A are essentially independent of the 
isospin asymmetry of the medium at high energies where the RIA is valid. Consequently, at these 
high energies the in-medium NN cross sections are also independent of the isospin asymmetry of 
the medium. Of course, there is still a difference between the neutron-proton and proton-proton 
(neutron-neutron) cross sections. In terms of these isospin independent quantities, the in-medium 
NN cross sections are obtained as 

a* nn = (A- 1 + X-^/pb, < p = (A- 1 - X-^/pb- (22) 

Here we assume that ct*„ = a* p , thus neglecting the small charge symmetry breaking effect [39] 
and the isospin-dependent Pauli blocking effects in asymmetric nuclear medium. The neglect of the 
former is also consistent with the model Lagrangian we used in this work. We do not consider the 
tiny proton-neutron mass splitting in free space. Moreover, our model Lagrangian similar to those 
used in most other RMF models does not include the scalar-isovector S meson responsible for the 
in-medium proton-neutron effective mass splitting. The same nn and pp scattering amplitudes are 
thus expected to result in the equal in-medium nn and pp cross sections. We stress here that the 
above results are only applicable to the high energy region where the RIA is valid and Pauli blocking 
effects are negligible. At low energies, the U sym and isospin-dependent Pauli blocking effects are not 
negligible, and the resulting NN cross sections will depend on the isospin asymmetry of the medium. 

III. RESULTS AND DISCUSSIONS 

The RMF model used here is the SL1 we constructed in our previous work [30]. The scalar 
densities used in the RIA are calculated with the SL1. The scalar densities calculated at low densities 
(p < L5po) arc quite similar with various models ranging from the best-fit models to the SL1. We 
refer the reader to our earlier work [30] for more details about the RMF with the SL1 parameter set. 

A. Schrodinger equivalent potential 

The proton and neutron SEP's obtained with the RIA serve as the basis to study the nucleon 
MFP and in-medium NN cross sections in this work. The isoscalar nucleon SEP, defined by the 
mean value of the neutron and proton SEP's in Eq. (21), is displayed in Fig. 1. It is seen that the 
imaginary SEP is stronger than the real one for nucleons at high kinetic energies. As the kinetic 
energy increases, the magnitude of the SEP increases very slowly. Consequently, this results in a 
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FIG. 1: (Color online) The isoscalar Schrodinger equivalent potential as a function of the nucleon kinetic 
energy for different densities. 

correspondingly small decrease of the nucleon MFP, as given by Eq.(17). Shown in Fig. 2 are the real 
and imaginary parts of the symmetry potential, defined by Eqs.(13) and (14), respectively. Indeed, 
we find that these potentials are independent of the isospin asymmetry of the medium. For the 
U sym , this was also found previously [33]. In Ref. [33], it was pointed out that the real part of 
the symmetry potential at fixed baryon densities is essentially a constant at high nucleon kinetic 
energies. As shown in the upper panel of Fig. 2, this is also true here. The U sym is negligibly 
small, actually very close to zero, in the whole energy range considered. This finding is consistent 
with the tendency of the energy dependence of the Lane potential extracted from nucleon-nucleus 
scatterings, which indicates that the Lane potential decreases with the increasing incident energy 
up to about 100 MeV above which there is still no data[20]. The smallness of U sym shown in the 
figure numerically justifies the approximation of using the kn instead of the k^ p in Eq.(20). The 
imaginary part of the symmetry potential, which differentiates proton and neutron absorptions in 
nuclear medium, displays a different dependence on the kinetic energy from that for the U sym . It 
is particularly interesting to note that its sign changes around Ekin ~ 630AfeV from positive to 
negative. This change can be reflected in the energy dependence of the nucleon MFP's. 
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FIG. 2: (Color online) The symmetry potential as a function of the nucleon kinetic energy for different 
densities. 

B. Nucleon mean free path 

The nucleon MFP's, calculated from Eq.(15) or (17), are shown in Fig. 3. We find that at higher 
kinetic energies the nucleon MFP changes little with the energy and it is much less sensitive to 
the isospin asymmetry. The rapid change of the nucleon MFP occurs around Ekin ~ 630MeV. 
It is interesting to see that the isospin-splitting between the neutron and proton MFP changes at 
Ekin ~ 630MeV. The nucleon MFP is determined by both the real and imaginary parts of the SEP 
as seen in Eq.(17). The difference between the real parts of the neutron and proton SEP's is tiny 
compared to the nucleon momentum at high kinetic energies. In fact, we see from Fig. 2 that the real 
part, i.e., U sym , is close to zero. Therefore, the change of the isospin-splitting between the neutron 
and proton MFP's can be attributed to the sign change of the imaginary part of the symmetry 
potential as shown in the lower panel of Fig. 2. This change for MFP's of neutrons and protons 
around 600 MeV/ A may have interesting experimental consequences and certainly deserves further 
studies. 

At Ekin > 700 MeV, the nucleon MFP is insensitive to the isospin asymmetry of the medium. This 
insensitivity is due to the fact that at high nucleon kinetic energies the W sym , responsible for the 
difference between the proton and neutron MFP's, is small compared to the W sep that dominates the 
contributions to the nucleon MFP's. At lower kinetic energies, however, appreciable sensitivity to the 
isospin asymmetry exists just because the imaginary part of the symmetry potential is comparably 
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FIG. 3: (Color online) The nucleon mean free path as a function of the nucleon kinetic energy for different 
densities. 

not so small. Comparing the two panels of Fig. 3, it is seen that the sensitivity of the nucleon MFP 
to the isospin asymmetry is reduced with the increasing density. This is also understandable because 
of the decreasing nucleon MFP's and the drop of the W sym to W sep ratio at higher densities. 

C. In-medium NN cross sections 

As discussed briefly in the text below Eq.(20), in the energy range where the RIA is valid, the in- 
medium NN cross sections <r* p and <7*„ pp are independent of the isospin asymmetry of the medium. 
They depend only on the density and the nucleon kinetic energy. Before presenting our results, it 
is worthwhile to add some more discussions. Usually, the in-medium NN cross section relies on the 
isospin asymmetry at low kinetic energies. Such dependence is introduced by the isospin dependent 
nucleon SEP through the real part of the momentum fc#, given by Eq.(16). In the non-relativistic 
limit, it has k R = (2ME kln - 2MU l sep )^ 2 with i = n,p. The fc£ and k p R can have considerable 
difference due to the non- negligible U symi which, for instance, can be found in Refs. [40, 41], at 
various isospin asymmetries. In fact, this is in also related to the Pauli blocking which contributes a 
p 2 / 3 -dependence in the kinetic part of the symmetry energy [21, 42-45]. At low energies the different 
Pauli blocking for neutrons and protons in asymmetric matter are important. They can essentially 
lead to the dependence of the in-medium NN cross sections on the isospin asymmetry of the medium. 
In the energy region where RIA is valid, the kinetic energies are much higher (roughly above 400MeV) 
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with the relative nucleon momentum fc 3> hp. The Pauli blocking thus becomes irrelevant. 




FIG. 4: (Color online) The in-medium NN cross sections as a function of the nucleon kinetic energy for 
different densities. The data for free NN cross sections are taken from Ref. [46]. 

The cr* p and <J* nnpp are calculated from the nucleon MFP's with Eqs.(19) and (22). As shown 
in Figs. 4 and 5, the in-medium NN cross sections increase with the density. At first glance, this 
feature seems to be in contrast to that occurring at lower energies where the Pauli blocking plays 
an important role in reducing the in-medium NN cross sections. However, at higher kinetic energies 
the in-medium NN cross sections are not necessarily constrained to be a descending function of 
density. In Fig. 4, as a check we see that the calculated NN cross sections at the very low density 
of O.lpo are very close to the free-space NN scattering data. Moreover, the in-medium NN cross 
sections are shown to be linear in density. This can be easily understood. As seen in Eqs.(20) and 
(22), the in- medium NN cross section are eventually determined by the SEP's which are quadratic 
in relativistic optical potentials, while the latter, given by Eq.(7), is linear in density. We note 
that the scalar density in Eq.(7) is linear in density in the relatively low density region where the 
RIA is applicable. Considering the density factor appearing in the denominator of Eq.(22), the in- 
medium NN cross section is therefore linearly proportional to the density. Interestingly, the ascending 
tendency of the in-medium cross section with density has also be seen at higher energies using some 
other approaches [7, 8, 10]. For instance, in Ref. [7], with the closed time path Green's function 
approach it was found that the cr*„ pp increases with density at E^in > 240 MeV. 

It is also seen in Figs. 4 and 5 that <7* p depends more strongly on the density than <r* n pp , while 
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FIG. 5: (Color online) The in-medium NN cross sections as a function of the density for various kinetic 
energies. 
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FIG. 6: (Color online) The ratio of a^ p to <Tnn,n P f° r various densities. 

the latter is more sensitive to the kinetic energy than the former. This different density and energy 
dependence for a* np and a* nn pp is actually a result of the distinction between the isovector (T=l) and 
isoscalar (T=0) interactions. Shown in Fig. 6 is the ratio cr^ p /a* n pp as a function of nucleon kinetic 
energy. The ratio is an ascending function of density. This feature is also different from that at 
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low and intermediate energies. It is clearly shown that the magnitude of cr* n pp exceeds <7* p around 
Ekin > 600MeV depending on the density. For a comparison, the experimental value in free space 
is also shown in Fig. 6. The experimental cr* p /cr* n pp ratio becomes smaller than 1 above about 
580 MeV. The general trend is qualitatively consistent with the calculations in the medium at the 
low density limit. It is worth noting that the in- medium NN cross sections obtained in the RIA rely 
little on the phenomenological freedoms except for those in the RMF model to calculate the scalar 
density. In fact, within the density region we considered in this study the scalar densities are almost 
the same for the SL1 model and the best-fit models, such as, the FSUGold [48]. 

IV. SUMMARY 

We have evaluated the relativistic optical potentials for neutrons and protons in asymmetric nuclear 
matter at high energies based on the relativistic impulse approximation with the empirical NN 
scattering amplitudes and the nuclear densities obtained in the relativistic mean field model. The 
Schrodinger equivalent potentials obtained from the relativistic optical potentials are then applied 
to calculate the nucleon mean free paths and the in-medium NN cross sections. The isospin-splitting 
of the proton and neutron mean free paths in asymmetric nuclear matter changes its sign at certain 
high energies, consistent with the changing sign of the underlying imaginary part of the symmetry 
potential. The small imaginary part of the symmetry potential at higher energies, compared to the 
counterpart of the isoscalar Schrodinger equivalent potential, accounts for the nucleon mean free 
paths with a much reduced sensitivity to the isospin asymmetry, while the negligibly small real part 
of the symmetry potential results in the in-medium NN cross sections independent of the isospin 
asymmetry of the medium. Moreover, in the energy range where the RIA is valid the in-medium NN 
cross sections are found to increase linearly with density. 
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